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ABSTRACT 

Objectives Lumbar disc degeneration (LDD) is prevalent, 
age-related and contributes to low back pain. Cross- 
sectional LDD as determined by MRI scan is known to be 
highly heritable. The authors postulated that the rate of 
progression might also be controlled by genetic factors. 
Methods A 10-year follow-up of MRI-determined LDD 
was performed in 234 pairs of twin volunteers in the 
UK and Australia, comprising 90 monozygotic pairs and 
1 44 dizygotic same-sex twin pairs. Of the total sample, 
95% were female. The mean age at baseline was 53.3 
years (range 32.3-69.5). The rate of progression was 
calculated and, because the effect of age was non-linear, 
the sample was divided into age strata and heritability 
estimated for each trait's progression. 
Results All MRI-determined traits worsened significantly 
over the period of follow-up (p<0.0001 for each). Change 
in disc height was not heritable at any age while posterior 
disc bulge was heritable across all age categories (range 
28-53%), with higher heritability in those over 60 years. 
Change in disc signal intensity and anterior osteophytes 
were found to be heritable only in those aged under 50 
years at baseline (heritability estimates 76% (95% CI 44% 
to 100%) and 74% (42% to 100%), respectively). 
Conclusions Longitudinal change in LDD traits is heritable 
for all traits except disc height, but there is a significant 
influence of age, which varies across traits. Future studies 
to define the genetic variants influencing LDD progression 
should examine MRI traits individually and in women should 
focus on those under 50 years of age. 



Low back pain is highly prevalent in the western 
world and accounts for considerable work absen- 
teeism. We and others have shown a relationship 
between back pain and lumbar disc degeneration 
(LDD), although the strength of association remains 
debated. 1-4 Several twin studies have shown LDD 
to be highly heritable but to date genetic analysis 
has focused on cross-sectional data. Previous lon- 
gitudinal analyses had suggested that progression 
may be more clinically relevant: the development 
of radiographic disc space narrowing over 9 years of 
follow-up, for example, was predicted by the pres- 
ence at baseline of low back pain in the Chingford 
study. 5 In this study, we examined LDD progres- 
sion longitudinally in monozygotic and dizygotic 
twins using lumbar spine MRI at two time points, 
over 10 years apart. We also estimated the herita- 
bility of LDD progression. 

SUBJECTS AND METHODS 

The study sample comprised same-sex twin pairs 
recruited from the TwinsUK register and the 
Sydney Twin Study. Between 1995 and 2001, 981 



male and female twins were subjects in the base- 
line MRI study. 6 None of the twins was aware 
of the hypothesis under investigation. The twins 
were neither selected nor excluded if there was a 
history of lower back pain, neck pain or degener- 
ate disc disease. This sample represents part of 
an adult female twin cohort recruited since 1995 
(www.twinsuk.ac.uk) from the general UK popu- 
lation 7 and from the Australian Twin Registry. 8 
The twins attended for MRI scan and a visit with 
a nurse metrologist to collect demographic data. 
Appropriate ethics permission had been obtained 
by both groups and twins gave informed written 
consent at both time points of the study. 

MRI at baseline 

The baseline MRI was performed using a Siemans 
LOT superconducting magnet (Siemans, Munich, 
Germany) in London and a GE Signa (GE, London, 
UK) 1.5T in Sydney. Sagittal images were obtained 
using a fast spin-echo sequence of repetition time/ 
echo time 5000-4500/112 ms with a slice thick- 
ness of 4 mm. 6 Serial sagittal images of the thora- 
columbar junction and lumbar spine (T9-L5) were 
obtained as described in detail, including reliability 
scores, elsewhere. At each disc level the main traits 
were scored by PNS considering the middle slices 
of the scan and using a progressive scale of 0-3 for: 
disc height measured in the middle of the disc; disc 
signal intensity within the nucleus pulposus; lum- 
bar disc extension posteriorly into the spinal canal 
and anterior osteophytes. A score of 0 indicates 
a normal disc, whereas a score of 3 is given for a 
highly degenerate disc. A summary score for the 
different traits was produced, by summing the cod- 
ing grades for the region Ll/2 to L5/S1. Twins of a 
pair were scanned on the same day, in the morning, 
more than 1 h after rising. 

MRI at follow-up 

The twins included in the baseline study were re- 
contacted between 2007 and 2009 and approxi- 
mately half the pairs re-attended for a clinical visit 
and further spine MRI scans. Scan parameters repro- 
duced those used at baseline (Figure 1). All imaging 
departments we contacted in London had upgraded 
their machines at some point over the interven- 
ing decade. Repeat scanning was performed on a 
Siemans 1.5T in London and a different Signa 1.5T 
in Sydney with the same slice measurements and 
excitation time and relaxation time as previously. 
Serial sagittal images of the thoracolumbar junction 
and lumbar spine (T9-L5) were obtained and coded 
for degenerative change using the atlas of Jarosz 
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et al, 9 as previously described in the baseline spine MRI study. 6 
MP was trained in coding by PNS who had coded the baseline 
scans. In brief, at each disc level the five middle slices of the sagit- 
tal MRI were used to score, on a 0-3 scale: disc height measured 
in the middle of the disc; disc signal intensity within the nucleus 
pulposus; lumbar disc extension into the spinal canal and anterior 
osteophytes. A score of 0 indicates a normal disc and a score of 
3 a highly degenerate disc. Scans were coded blindly, without 
reference to baseline coding, scan or to co-twin's results. For each 
trait, the baseline score was subtracted from the follow-up score 
to give a change in score. Because the interval between the two 
scans varied among the twins, the score was adjusted for the 
time interval between the two MRI scans — effectively giving a 
rate of change per year. Reproducibility scores were performed 
intracoder (mean of the four traits 0.68; comprising disc height 
0.77, disc signal 0.61, disc bulge 0.67, anterior osteophytes 0.66) 
and intercoder (mean 0.54; disc height 0.66, disc signal 0.62, disc 
bulge 0.56, anterior osteophytes 0.32). 

Statistical methods 

Comparisons were made using t tests, the Wilcoxon signed rank 
test and multivariate regression analysis as appropriate using 
Stata 2007. Intraclass correlations were calculated for monozy- 
gotic and dizygotic twins. The underlying distribution of change 
in score is considered to reflect contributions by additive genetic 
(A) factors, shared environmental factors (C) and additive envi- 
ronmental factors (E), and the influence of each was modelled 
using the De Fries Fulker regression approach implemented in 
Stata. The effect of age on the rate of progression was non-linear 
in this dataset so the analysis was carried out in three baseline 
age strata: 50 years or less; 50-60 years; 60 years or over. 

RESULTS 

Longitudinal data were available for 234 same-sex twin pairs 
of confirmed zygosity, comprising 90 pairs of monozygotic and 



144 pairs of dizygotic twins. The majority (211 pairs) were from 
the TwinsUK cohort, the remaining 23 pairs were recruited in 
Australia. Twenty-four subjects (5.1%) were male, distributed 
evenly between the two zygosities. Table 1 shows the charac- 
teristics of the total sample by zygosity. Mean age at baseline 
was 53.6 years (range 40.1-68.7), with mean age in monozy- 
gotic twins 55.4 years and dizygotic twins 52.6 years. The 
mean interval between baseline and follow-up was 10.7 years 
(range 7.6-13.7). The Australian twins were significantly older 
(Australian vs UK 56.2 vs 53.4 years, p=0.01) and heavier (body 
mass index (BMI) 27.6 vs 24.5 kg/m 2 , p<0.01) than the UK twin 
volunteers, but there was no difference in baseline levels of LDD 
between the two groups. In general, the monozygotic twins had 
more advanced degenerative change at baseline (table 1), with 
a significantly greater osteophyte score and disc signal in the 
monozygotic compared with the dizygotic twins. These differ- 
ences were attributable to age: when considered in age strata the 
baseline traits did not differ significantly between the zygosi- 
ties (table 1). There was no significant difference between twins 
attending for follow-up and those not attending in terms of age 
and baseline summary LDD score but there was a difference 
in one of the MRI subtraits, osteophytes, with those not fol- 
lowed up having significantly greater baseline levels than those 
re-attending (1.79 vs 1.14, respectively, p<0.001). 

Progression of LDD 

All MRI traits showed an increase in score over the period of 
follow-up (table 2), indicative of a deterioration in LDD. In 
fact, almost all individuals demonstrated deterioration over the 
period of follow-up: only one person (0.2%) had the same sum- 
mary LDD score at both time points and one person (0.2 %) had 
a score that had marginally improved. The change in summary 
score over the period of follow-up was highly significant, the 
mean LDD summary score at baseline versus follow-up 12.96 
versus 22.97 (p<0.0001). Dizygotic twins were found to have 



Table 1 Sample characteristics at baseline by zygosity, showing mean MRI trait scores for whole group and 
by age strata 





Age, years 


Monozygotic twins 


Dizygotic twins 


p Value 


Sample size 




180 


288 




Age, years, mean (range) 




55.42 (40.89-68.74) 


52.64 (40.08-68.75) 


<0.01 


Males, n (%) 




12 


12 


0.233 


Body mass index, kg/m 2 (range) 




24.39(18.85-42.31) 


25.01 (17.51-45.14) 


0.10 


Lumbar disc height score, mean (range) 


All ages 


4.07(0-13) 


3.80(0-12) 


0.25 




<50 


3.08(1.59) 


3.48 (2.24) 


0.33 




50-60 


3.95 (2.44) 


3.80(3.36) 


0.64 




>60 


5.2142 (3.15) 


4.49(2.71) 


0.23 


Lumbar disc signal score, mean (range) 


All ages 


6.21 (0-11) 


5.65(0-12) 


0.03 




<50 


4.19(2.29) 


4.64(2.69) 


0.37 




50-60 


6.37 (2.44) 


5.80(2.60) 


0.09 




>60 


7.57 (2.20) 


7.41 (2.15) 


0.73 


Lumbar disc bulge score, mean (range) 


All ages 


2.15(0-8) 


1.97(0-10) 


0.24 




<50 


1.83(1.42) 


1.87(1.64) 


0.89 




50-60 


2.12(1.707) 


1.84(1.54) 


0.20 




>60 


2.5(1.81) 


2.47(1.77) 


0.94 


Lumbar anterior osteophyte score, mean (range) 


All ages 


1 .40 (0-7) 


0.99(0-10) 


<0.01 




<50 


1.14(1.59) 


0.77(1.11) 


0.13 




50-60 


1.24(1.16) 


0.96 (1.22) 


0.08 




>60 


2.0(1.82) 


1.53(2.06) 


0.25 


Summary lumbar score, mean (range) 


All ages 


13.84(0-38) 


12.41 (0-40) 


0.02 




<50 


10.25(4.93) 


10.77(5.78) 


0.63 




50-60 


13.68(6.08) 


12.39(6.05) 


0.12 




>60 


17.29 (6.55) 


15.91 (6.86) 


0.32 



Summary lumbar score represente for all four traits over five lumbar disc levels, p Value is significance value comparing monozygotic 
with dizygotic twins. 
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Table 2 Changes in LDD phenotypes over time, by zygosity 





Monozygotic twins 




Dizygotic twins 






Comparing monozygotic and dizygotic 




Baseline 


Follow-up 


Change in 


Baseline 


Follow-up 


Change in 


p Value 


p Value 


p Value 


Trait 


score (SD) 


score (SD) 


mean score (SD) 


score 


score 


mean score (SD) 


baseline 


follow-up 


change in trait 


Disc height 


4.07 (2.58) 


6.39 (2.65) 


2.27 1.93 


3.80(2.44) 


6.52 (2.82) 


2.71 (2.14) 


0.25 


0.63 


0.03 


Disc signal 


6.21 (2.61| 


8.55(2.13) 


2.31 2.25 


5.65(2.73) 


8.25 (1.98) 


2.62(2.18) 


0.03 


0.13 


0.15 


Disc bulge 


2.15(1.69) 


3.88(2.37) 


1.65 (1.76) 


1.97 (1.63) 


4.01 (2.43) 


2.07 (2.00) 


0.24 


0.55 


0.02 


Anterior osteophytes 


1.40(1.46) 


4.22 (2.68) 


2.72 2.40 


0.99(1.40) 


4.15(2.70) 


3.16(2.47) 


<0.01 


0.76 


0.06 


LDD summary score 


13.84 (6.39) 


23.05 (6.88) 


8.95 4.27 


12.41 (6.35) 


22.93 (7.53) 


10.55* (4.68) 


0.02 


0.86 


<0.01 



Mean summary scores for the individual traits over the five lumbar discs and SD are shown for baseline and follow-up scans, with mean change in score over the 1 0-year period of 
follow-up, and p values are shown for comparison between monozygotic and dizygotic twins at baseline, follow-up and for change in trait. Taken together, all traits' changes over 
10 years were highly significant, each trait p<0.0001. 
LDD, lumbar disc degeneration. 



Table 3 


Intraclass correlations, modelling and heritability estimates for the MRI traits by age category 






Age, 




N pairs 


N pairs 


R-monozygotic 


R-dizygotic 






years 


ATrait/time (per year) 


(monozygotic) 


(dizygotic) 


(95% CI) 


(95% CI) 


A (95% CI) 


C (95% CI) 


<50 


Disc height 


18 


53 


0.51 (0.16 to 0.87) 


0.53 (0.33 to 0.72) 




0.51 (0.30 to 0.71) 


50-60 


Disc height 


50 


63 


0.33 (0.076 to 0.57) 


0.40 (0.1 9 to 0.62) 




0.37 (0.19 to 0.54) 


>60 


Disc height 


22 


28 


0.18 (0.00 to 0.63) 


0.47 (0.17 to 0.77) 




0.38 (0.10 to 0.65) 


<50 


Disc signal 


18 


53 


0.71 (0.48 to 0.95) 


0.41 (0.18 to 0.63) 


0.76(0.44 to 1.0) 




50-60 


Disc signal 


50 


63 


0.37 (0.13 to 0.61) 


0.34(0.11 to 0.57) 




0.36 (0.18 to 0.54) 


>60 


Disc signal 


22 


28 


0.45 (0.08 to 0.81) 


0.48 (0.18 to 0.78) 




0.41 (0.14 to 0.68) 


<50 


Disc bulge 


18 


53 


0.39 (0.0 to 0.80) 


0.10 (0.00 to 0.37) 


0.28(0 to 0.65) 




50-60 


Disc bulge 


50 


63 


0.56 (0.38 to 0.75) 


0.16 (0.00 to 0.41) 


0.51 (0.28 to 0.74) 




>60 


Disc bulge 


22 


28 


0.55 (0.24 to 0.86) 


0.34(0.00 to 0.68) 


0.53 (0.13 to 0.94) 




<50 


Anterior osteophytes 


18 


53 


0.71 (0.49 to 0.94) 


0.40 (0.1 8 to 0.63) 


0.74(0.42 to 1.0) 




50-60 


Anterior osteophytes 


50 


63 


0.58 (0.40 to 0.77) 


0.58 (0.42 to 0.75) 




0.58(0.43 to 0.74) 


>60 


Anterior osteophytes 


22 


28 


0.60 (0.31 to 0.89) 


0.56 (0.29 to 0.82) 




0.57(0.32 to 0.81) 


<50 


Summary score 


18 


53 


0.57 (0.24 to 0.89) 


0.12 (0.00 to 0.39) 


0.40 (0.06 to 0.73) 




50-60 


Summary score 


50 


63 


0.42 (0.19 to 0.65) 


0.31 (0.08 to 0.55) 


0.48 (0.24 to 0.72) 




>60 


Summary score 


22 


28 


0.39 (0.00 to 0.77) 


0.48 (0.1 8 to 0.78) 




0.40 (0.13 to 0.67) 



Heritability was calculated using 90 monozygotic and 1 44 dizygotic same-sex twin pairs (omitted when dizygotic intraclass correlation greater than monozygotic intraclass correlation 
clearly suggests no heritable component) when N is the sample size in twin pairs by age subgroup. Age category 0=age <50 years, category 1 =age. ATrait/time per year represents 
change in trait per year. A represents additive genetic, C common environment and E unique environmental factors. 



significantly greater deterioration in LDD than monozygotic 
twins (10.55 vs 8.95, p=0.0002) and, adjusting for the time 
interval between the two MRI scans, a greater rate of progres- 
sion (1.0 vs 0.84 per year, p<0.001). BMI, age at baseline and 
gender had no detectable effect on the progression of LDD over 
10.7 years. 

Heritability of progression of LDD 

Table 3 shows the intraclass correlation coefficients for the 
monozygotic and dizygotic twins by MRI trait, together with 
the associated heritability estimates derived from modelling. 
Intraclass correlations being greater in monozygotic than dizy- 
gotic twins is indicative of a genetic influence on a trait. Disc 
height progression showed no evidence of influence by genetic 
factors in any age stratum. In the other traits there was evidence 
of genetic influence particularly in the younger age group: disc 
signal in those less than 50 years had a high heritability esti- 
mate of 76%. In contrast, disc bulge was found to be heritable 
at all three age categories, with greater heritability in older sub- 
jects, range 28 (0-0.65%) to 53 (0.13-0.94%). The progression 
of osteophytes in those less than 50 years showed a marked 
genetic effect (0.74 (0.42-1.05)) with no genetic influence detect- 
able in the older age categories. Unsurprisingly, the heritability 
of the summary score progression reflects that of its components 
and is mainly influenced by genetic factors at younger age. 

DISCUSSION 

This is the first report of the heritability of progression of LDD 
and we have an extended period of follow-up of over 10 years 



in our sample of monozygotic and dizygotic twins. The main 
findings were that all the MRI-determined traits deteriorated 
over the decade and genetic influences on progression varied 
across the traits as well as by age. For example, change in disc 
height did not appear to be influenced by genetic factors at 
any age while disc bulge was heritable in all age categories. 
Deterioration of disc signal intensity and osteophytes had a 
significant genetic influence at young age only, which may be 
a reflection of the greater variation in progression in this age 
group. We previously reported no significant heritability for 
disc signal intensity based on the cross-sectional (baseline) MRI 
findings in this cohort. 6 Taken together, these data suggest that 
any genetic effect on disc signal occurs early rather than late in 
adult life. This is not without precedent: heritability of a num- 
ber of other traits, including musculoskeletal traits, is known 
to vary with age. For example, there is no detectable genetic 
contribution to hip fracture in older people (over 80 years), 
whereas the trait is clearly heritable in younger people, with 
estimates from a large twin study in Sweden suggesting herita- 
bility as high as 68%. 10 

The progression of LDD over a 5-year time period has been 
investigated previously by plain radiography in UK singletons, 5 
which showed a 4-5% progression rate, and by MRI in 134 male 
monozygotic twins from Finland showing 3-4% progression. 11 
The latter study highlighted different MRI-determined traits 
progressing at different rates. As the Finnish study did not 
include dizygotic twins, the heritability of progression could not 
be estimated but the investigators used the less specific measure 
of 'familial aggregation' determined from intraclass correlations. 
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Twin l Twin 2 



Baseline 



Follow-up 




Figure 1 Illustrative MRI lumbar spine scans from a pair of 
monozygotic twins. Sagittal T2-weighted MRI scans from a pair of 
female monozygotic twins, baseline scan uppermost. Aged 52 years 
at baseline, both follow-up scans show increasing loss of disc signal 
intensity, loss of disc height and increasing extent and severity of disc 
bulges posteriorly. 

Their findings of MRI-determined posterior disc bulge progres- 
sion not being genetic are different to ours, as were their results 
for disc height. 11 However, the Finnish study was of a consid- 
erably smaller sample (n=67 twin pairs; 3.5 times smaller than 
the sample reported here) over a shorter time period (5 years) 
and the authors did not study changes in signal intensity. The 
other major difference between the two cohorts is in the gender 
of the participants, so perhaps the contrasting findings reveal a 
true sex difference in LDD progression. Our findings highlight 
the heterogeneity of the MRI traits, and suggest that they reflect 
very different features of the pathological process and are under 
differential control. This suggests that future studies should con- 
tinue to examine these features separately. This is in keeping 
with findings from the osteoarthritis literature — the approach 
is commonly referred to as a 'splitting' rather than 'lumping' 
method. These terms have also been used to describe the practice 
of summing osteoarthritis across body sites rather than consid- 
ering each individually. Support for this approach is mounting, 
with the accumulating evidence of differential genetic influences 
at different anatomical sites. 12 13 

We acknowledge a number of limitations to this study. As 
our sample was predominantly of female subjects it does not, 
on the whole, contain subjects reporting heavy manual occu- 
pations, although such environmental factors are now consid- 
ered to play a small role. LDD is thought to be similar across 
the sexes; although site and sex-specific genetic predisposi- 
tion is being identified for peripheral joint osteoarthritis, there 
is little evidence of a significant structural difference between 
the sexes in LDD. Men are thought to have earlier degenerative 
change based on a postmortem study from the 1920s, 14 but as 
we detected no significant difference in LDD between the men 



and the women in our sample, we felt it reasonable to include 
both groups in our study. These twins have been shown to be 
representative of the UK population for a range of traits includ- 
ing degenerative joint diseases. 15 Twins were not recruited on 
the basis of back pain but the possibility that this may have 
been a motivation to participate cannot be fully discounted. A 
concern in these data was that there was a relatively poor mea- 
sure of intercoder reproducibility (k=0.54) and systematic dif- 
ferences in the monozygotic and dizygotic groups. The latter is 
partly the result of slightly younger dizygotic twin volunteers, 
and while this is probably a chance finding, it is possible that a 
degree of selection bias introduced systematic differences. The 
intercoder K is not likely to have influenced findings as almost all 
the baseline scans were coded by PNS and the follow-up scans 
by MP. Systematic bias in the scans and their coding could have 
influenced the rate of progression of LDD — newer MRI scans 
were inevitably of higher quality than the baseline scan — but 
one would not expect a differential between monozygotic and 
dizygotic twins, so heritability would not be affected. 

Our initial data analysis showed that the effect of age on the 
rate of degeneration was non-linear. Furthermore, there was 
a suggestion of a ceiling effect in the coding for degenerative 
change, with more severe disease progressing less rapidly. As a 
result, we felt that it was most appropriate to examine herita- 
bility in age strata, rather than by attempting to eliminate age 
effects using regression modelling. This approach necessarily 
reduced the statistical power, and is reflected in the wide CI of 
the stratum-specific heritability estimates (table 3). While defini- 
tive statements regarding the genetic influence on the progres- 
sion of MRI traits at differing age groups should be made with 
caution, the results across traits are suggestive of a change in 
genetic influence over age. 

That genetic factors underlie the rate of progression of LDD is 
not, perhaps, surprising in itself, but is of importance because it 
provides the rationale to the search for precise genetic variants 
controlling LDD progression. In addition, progression of LDD 
may be more closely related to symptoms of low back pain than 
cross-sectional degeneration — there is evidence that this is the 
case in similar phenotypes such as knee osteoarthritis. 16 These 
results will direct the search for genetic variants influencing the 
rate of deterioration as, arguably, those at risk of rapid progres- 
sion are the most important group to identify when interven- 
tions become available for trials. Cross-sectional LDD is highly 
heritable and has been the subject of a number of candidate gene 
studies and variants identified that have been replicated in one 
or more studies include VDR, COL9A2, COL9A3, COL1A1, 
MATN3, MMP and some inflammatory genes such as IL-1 and 
THSD2. 17 However, in common with many complex traits, the 
genetic variants described thus far account for only a small pro- 
portion of the estimated heritability, and it is likely that many 
will turn out to be false positives, as has been reported for 
osteoporosis. 18 

Little is known of the environmental factors governing the 
rate of progression of LDD. BMI is strongly suspected if not 
proved, smoking and occupational maximum lifting have been 
shown to influence progression in twin studies. 19 Good longitu- 
dinal epidemiological studies remain scarce (for cross-sectional 
review see Williams and Sambrook). 20 This information will 
be important information to obtain, for if the main effects are 
environmental then these may be modifiable by, for example, 
weight loss or occupational changes. If, however, most of the 
variation in progression is attributable to genetic factors, then 
further understanding of gene products may provide novel 
targets for intervention. In addition, the identification of new 
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biological pathways may aid in the identification of much 
needed biomarkers for future trials of disease-modifying drugs. 
Low back pain is a considerable public health issue: new path- 
ways to target would be of great interest to the pharmaceutical 
industry and to the many practitioners with a role in managing 
back pain. Our results point the way to subject selection in stud- 
ies to identify genetic variants: they highlight that such changes 
are largely seen in young-middle age. Genetic variants influenc- 
ing progression are now a realistic target for researchers into 
LDD. This study provides the impetus for further work to iden- 
tify the genetic variants involved as well as a clinical assessment 
of whether those with rapid LDD progression are at greater risk 
of symptoms. 
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